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Abstract. Incompressible fields are of special importance in elec- 
trodynamics, fluid mechanics, and quantum mechanics. We shall derive a few 
| expressions for such fields in a Riemannian manifold, and show how to generate 

an incompressible field from an arbitrary set of scalar differentiable functions. 
The concept of compressibility removing factors of an arbitrary vector field 
is introduced and utilized to obtain from an arbitrary vector field an incom- 
| pressible one that has the same vector surfaces as the original field. A general 

expression for compressibility removing factors of a vector field is derived. The 
. method is applied to central fields. 

o . 

^ ■ 1. Introduction 

The equation which expresses that a vector field L has a zero divergence, divL = 0, 
o r appears in a number of areas of physics such as electrodynamics, fluid mechanics, and 

quantum mechanics. In fact two of Maxwell's equations, namely div E— p/to and 

. -> — » — » 

div H= 0, where E and H are the electric and magnetic fields respectively, p is the 
charge density and eo is the permittivity of the medium, reduce in the free space to 

the form we have mentioned. Also the continuity equation div j +dp/dt = 0] in 

which j is the current density of some flow of matter or charge, p is its volume density 
and t is the time, is reducible to this form when the volume density does not change 
with time. More recently ||, it was shown that the incompressible vector fields in 
a manifold, single out a class of quantum momentum observables that exhibit the 
features of the Cartesian momentum operators in a Euclidean space. In fact it was 
shown that the eigenfunctions of any member of this class is a direct generalization 
of De Broglie waves. The solution of the equation divL = in the three dimensional 
Euclidean space is well known. It is given in terms of a C 2 vector field A, by L= 
Vx A ■ In this work 

(1) We shall find a few expressions for the solutions of this equation in an n- 
dimensional Riemannian manifold. 

(2) Show that every ,(n — 1),C 2 functions determine an incompressible field. 
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(3) Introduce the concept of a compressibility removing factor of a vector field, 
construct some methods for determining such factors and find the most general ex- 
pression of such factors. 

(4) Show that the quotient of two compressibility removing factors of a field is a 
vector surface of that field. 

(5) Apply the results we have obtained to central fields and centrally symmetric 
fields. 

We shall use tensor notations in aim to simplify the calculations and formulae 
obtained. However the reader must be warned that the expressions we get for in- 
compressible vector fields, though valid in any coordinate system, have no tensorial 
character, i.e. they do not transform as a contravariant vector. 

2. The Field Equation 

Let M n be a Riemannian manifold with ( X ■ .... j OC ^ 3jS 3j chart in which the metric 
takes the form ds 2 = gijdx l dx^ , and let g = det(gijy Let L be a C 1 vector field 
on M n which has in the given chart the form L = C, l d/dx l .The field L is said to be 
incompressible if its divergence 

divL = g- 1 l\g 1 ' 2 C), (1) 
is zero [[jj]. According to (P the field L is incompressible if 

(V9?),i = 0- (2) 
The incompressible field equation (0) may be written as 

A\ = o where A 1 = y/gC- 

Equation (||) admits obvious solutions such as 

C = g~ 1/2 , ? ()(./ 1 //:.///•). (3) 

where r is a fixed index chosen arbitrarily from the set {1,2, ....,n}. Since the field 
equation (§) is linear in the unknown functions ^(orA^) all linear combinations in 
the solutions (ff) are also solutions of the field equation. We shall see also that 
the commutator of two incompressible fields is also incompressible, i.e., the set of 
incompressible fields in M n forms a Lie algebra. 

In the case of an one dimensional manifold M\ which is coordinated by x so 
that ds 2 = g(x)dx 2 , the most general solution of (0) is L — -j=-^ , where c is an 
arbitrary constant. However Mi is certainly Euclidean and the metric in Mi may be 

reduced to the standard form ds 2 = dX 2 , by the transformation X =J ^fgdx. The 
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incompressible field in the new coordinate X takes the form L = c d/dX which is 
an infinitesimal motion of M%. When M\ is isometric to the real line $ft or to the one 
dimensional circle S , the field L generates the symmetry group U(X) — X + t, (t G 
of motions of Mi. 

In a two-dimensional manifold M2, a field L is incompressible if A\ + A 2 2 = 0, a 
condition which is satisfied if and only if the ordinary differential equation A 2 dx 1 — 
A 1 dx 2 = is exact; which in turn is equivalent to the existence of a function B of 
class C 2 on M 2 such that B \ = —A 2 and B 2 — A\. It follows that the most general 
form of an incompressible field in M 2 is of the form 

1 dB d d# (A) 

y/g dx 2 dx 1 dx 1 dx 2 

It is tempting to generalize the last result in a straight forward manner to an M n . 
To achieve that we start with an arbitrary C n function B, and take A = y/gC given 
by 

A = C\B 2Z...n > A = C2-Bl3.... n , , A™ = C n B ,l2....n-l) (5) 

n 

where c, are arbitrary constants such that q = 0. We have then 

i=i 

n 

4 =E ^ 5 ,12....n = 0, 

i=l 

which proves that the field L, as given by (H), is incompressible. If all components of 
L are non-zero we may satisfy Q — by 



-il- 



where no summation over i is implied in the last equation. When the dimension 
of the manifold is even one may choose an alternative solution c, = (— This 
method however, determines only a narrow class of incompressible fields in M n . i.e., 
if L — ^d/dx 1 is incompressible then there may not exist a function B such that @ 
is satisfied. The following theorem singles out the case in which an incompressible 
field L may be expressed by (0) 

Theorem 1. An incompressible field L = t^d/dx 1 can be expressed by the form ([|) 
if and only if there exist constants Cj(i = 1, n) such that c i — and 

-A\ = -A% = .... = -A%. (6) 

Ci ' C 2 ' Cn ' 

If for some index r, A r r = then we take c r = and exclude the corresponding ratio 
from the equalities (&). 



For proof see appendix 1. 
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3. General Expressions for an Incompressible Field 

3.1. A First Expression:. This results from taking {n — 1) components of L 
as arbitrary functions of class C l and calculating the remaining component £ r (r ^ j) 
by the field equation (|2]); it is given by 

v^f = "E J(V9e),dx r + h r , (7) 

where h r is an arbitrary function of class C l independent of x r . Choosing every 
component £ J as the partial derivative with respect to x r of an arbitrary function 
of class C 2 , we write ) in the equivalent form 

y/g& = 4 U * r) , y/g? = - E s j- 

We note that the arbitrary function h r has been absorbed in B\. 

3.2. A second Expression:. Let e* 1 be the totally antisymmetric unit tensor 
I 

e* 1 = 1 if (ii....i n ) is an even permutation of (1,2, ...,n) 
= — 1 if (ii....i n ) is an odd permutation of (1, 2, n) 
= if two indices are equal, 

and let B h ,,, An _ 2 , where i n _ 2 G {l,...,n} and i± < i 2 < .... < z n -2 , be 

arbitrary (™) functions of class C 2 on M n . 

Theorem 2. Tie field 

L = -U' 1 i -B il ... An _ 2 , in _ 1 d/dx i - (8) 

is incompressible. 
Proof: See appendix A2 

3.3. A Third Expression:. Let 

<t> k {x\....,x n ) fce[l,n-l] (9) 
be (n — 1) arbitrary C 2 functions on M n that are functionally independent, i.e. 

rank(d(p k /dx^) = n — 1 
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Theorem 3. The field 



L = ^ ix - An Mn, <t>n-i,i-i-^~ (10) 



in M n is incompressible. 
Proof: See appendix A3. 



Remark 1. The expression ( [it]) may be written as 



1 

71 



'1,2 



'1,71 



<Pn-l,l </>n-l,2 <Pn-l,n 

d/dx 1 d/dx 2 d/dx n 



(11) 



where the determinant must be expanded in a way that the minor of d/dx 1 is the 
component of £\ 

Remark 2. Had the functions (|^) were functionally dependent then 

rank(d(pk/dx' 1 ) < n — 1, k G [1, n — 1], «6 [l,n], 

and the determinant of every (n — 1) x {n — 1) sub-matrix of the last matrix vanishes. 
However these determinants, apart from a possible difference in sign, are precisely 
the components of L. Hence L vanishes. 



Theorem 4. Each of the surfaces 

cf )k (x\....,x n ) = c k fc 6 [1,71-1] (12) 

is a vector surface of the field ( fltl ), and hence the intersection of these surfaces are 
the field's lines. 

Proof. Let the field L, as given by the expression (|iT|), acts on It is apparent 
that the determinant on the right hand-side vanishes, on the account of its k-th and 
last rows are equal. 

Remarks and Examples 

1. We have seen that in M 2 there exists a unique expression of an incompressible 
field, namely (£|). Applying any of the last expressions we get @) 
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2. When we apply the second expression to M 3 we get 

L = g- 1/2 e ilk B hJ d/dx k 

= <T 1/2 [(£ 2 , 3 - B^d/dx 1 + (£?3,i - B^)d/dx 2 + (S li2 - 5 2 ,i)<9/«9x 3 ] 

where Bi,B 2 and S3 are arbitrary functions of class C 2 . This expression is equiva- 
lent to the familiar expression L= Vx B of an incompressible field in a rectangular 
coordinates in a Euclidean manifold M3. 

Applying the third expression in M3 we get 



9 1 V- M o,r„ = ^-V 2 



</>,l </>,2 </>,3 

^,2 ^,3 

d/dx 1 d/dx 2 d/dx 3 



This reduces when M 3 is Euclidean and the coordinates are Cartesian to L= V0x V?/>. 

3. The expression (|l(]) of an incompressible field L is equivalent to obtaining the 
system of ordinary differential equations, whose integrals , k G [1, n — 1] are given, 
through taking the differentials d(f>k = 0, solving for (n — 1) of the dx l ,a.nd eventually 
symmetrizing the system of ordinary differential equation that has been obtained. 

4. Reduction of a Vector Field to an Incompressible One 
Let K = rfdjdx 1 be a C 1 vector field on M n , and n(x l , ...,x n ) be a C 1 function on 
M n . The function /1 is said to be a compressibility removing factor (CRF for short) 
of the vector field K, or it removes the compressibility of K, if the vector field \xK 
is incompressible. If /i removes the compressibility of K we say that fiK reduces 
the field K. A function /i is a CRF of K if and only if it satisfies the equation 
div(fiK) = 0, which is equivalent to fidivK + K\i = 0. This equation can be deduced 
directly, or by using the properties of div [|I|. The equation of CRF 

Kfi = —fidivK (13) 

is a partial differential equation of Lagrange type. The subsidiary system associated 
with this first order linear equation [|5| is 

dx 1 dx 11 da , . 

= = = (14) 

r/ 1 r/ n — /i divK 

We shall denote the set of incompressible fields which reduces K by L(/i, K). 



Theorem 5. Each member of the set L(fi, K) has the same vector surfaces of K. 
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Proof. It follows from observing that the vector surfaces of all the fields u,K, 
as well as the field K, are determined by the first n ratios in (|14]) . Or from noting 
that fjJC<f> = <£> K(p = 0. 

Remark 3. The last theorem is applicable to every vector field of the form XK, 
where A is any function that is not zero. This is true whether or not A reduces K. 

Let 

k (x\....,x n ) =c k fee [1,71-1] (15) 



be functionally independent integrals of ([14]) resulting from the first n ratios. The 
most general arbitrary integral associated with the field K is of the form f((j)\, n -i) ; 
where / is an arbitrary function of the Hence every vector surface of the field K 
has the form f(<fri, 4> n -i) = 

Theorem 6. If the CRFs \L\ and \x 2 of the field K are linearly independent then 
/ii//i 2 is an integral of K. 



Proof. By the equation of CRFs (|T3"D, ^K^i = ^\K^2 , and hence 

K(jJ,i/fi2) = /i 2 2 (fi 2 Kfi 1 - fiiKfi 2 ) = 

Therefore ^1/^2 is an integral of the subsidiary system associated with K, and 
^1/^2 = c, where c is an arbitrary constant, is a vector surface of K. 

Theorem 7. If fi is a CRF of the field K then all its CRFs are given by 

M = /i/(0 1 ,....,0 n _ 1 ), (16) 

where f is a C l arbitrary function in the functions 

Proof. Let M be a CRF of the field K. By theorem (4.2) M/fi is an integral 
of the subsidiary system associated with the field K, and hence it has the form 

M/fi = /(01,....,0 n _l), 

which completes the proof. 

One can verify easily that the equation of CRFs ( |i~3|) is satisfied by M given by 

(0): 



(nf)divK + K(/if) = fifdivK + fKfj, + \iKf 

= f.QidivK + Kfjt) + (becase Kf = 0). 
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By ([13]) the right hand side vanishes. 

It follows from theorem (^j that if K is incompressible then every CRF of K is 
of the form M = f(<pi, <f> n -i)- 

We have seen in section 2 that a field L = ^d/dx 1 + £ 2 d/dx 2 is incompressible if 
and only if the equation 

v/^W - ^dx 2 = (17) 

is exact. It follows that \x is a CRF of L if and only if \x is an integrating factor of 
(|T7|), i.e. if and only if /x-y/g is an integrating factor of £ 2 dx l — ^dx 2 = o. 

We finally mention that, although there exists an infinite family L(/x, K) of fields 
which have the same vector surfaces as K does, the integral curves of these fields are 
distinguishable. Indeed, the differential equations of the integral curves of a field \iK, 
namely 

dx l = fit] 1 dt, 

depend on the function /x. 

5. Method for Determining CRFs 
1. If it were possible to find an integral of the system (|14|) that contains /x then on 
solving for /x we get one of the required factors. 

We shall assume in what follows that we have available (n — 1) integrals of the 
field K which we denote by <pk(k = 1, n — 1). 

2. Since the integrals <pk of K are also integrals to every incompressible field 
L = fiK, we may find one of the fields L = fiK say L = t^d/dx 1 using the expression 
(|lOf) of an incompressible field in terms of its vector surfaces. Assuming that the 
component rf of the field K is not zero, then 

H = C/ V r = e il -^- ir 0i )il ....0n-i )in _ 1 A7 r (18) 

is one of the required factors. 

3. Assuming that rank(d<f)j/dx k ) = n — 1 where k,j G [i, n — 1], we may solve 
(p~5| ) for (x , ....,x n ~ l ) in the form 

x k = x h (c 1 , ...,c n _i,x n ) fce[l,n — 1]. (19) 

If we could obtain the solution (|T9"|) we substitute in (|T^) , separate the variables and 
integrate to find 

,x = Cn e-f^ dxn = c n e- F ^--> c "-^ n \ (20) 

where x k , k e [1, n — 1], are substituted for in F = / ^rdx n from (JT5|) . Substituting 
for Ci, c n _i in ( ^0|) from (|19"D we get the n-th integral of (|T^) 



/ x(x 1 ,....,x") = /(0 1 ,...,0 n _ 1 )e- F ^' >^ n \ 
which is also the general form of a CRF. 
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6. Example - Central Fields 
Let M n = E 3 be the three dimensional Euclidean space with Cartesian coordinates 
r= (x 1 = x, x 2 = y,x 3 = z). We shall say that the field K = r\ l djdx % is central if 
its lines are straight lines through the origin; and centrally symmetric if it can be 
reduced to the form K = rj(r)d/dr, where r = r . 
Consider the centrally symmetric field 

K = xd/dx + yd/dy + zd/dz. (21) 

Since divK = 3, the equation of CRFs is 

Kn + 3// = 0. (22) 

The subsidiary system associated with ( ^2j) is 

dx/x = dy/y = dz/z = —dfi/3fi. (23) 



From (|23"D we deduce that each of the functions /ii = x~ 3 ,/i2 = y~ 3 ,fis = z~ 3 is a 
CRF of K. It follows that /Ui//i 2 = y 3 /x 3 , /ii//i3 = z 3 /x 3 are two integrals of the 
subsidiary system associated with K, and hence every vector surface of K is of the 
form 

f(y/x,z/x) = 0. (24) 
According to (JTB|) the general form of a CRF is 



fi = x 3 f(y/x,z/x). (25) 

We may choose / so that we get the CRFs /i 4 = 1/xyz or fi 5 = (x 2 + y 2 + z 2 ) 3 ^ 2 . 
Every vector surface of the fields \iiK{% = 1,...,5) is given by (|2HD, which shows 
that the lines of all these fields are straight lines through the origin, and hence 
these fields are central. The field K in spherical coordinates (r, 0, 9) is written as 
K = rd/dr which shows that it is centrally symmetric, and that the product of K by 
any /ij, except /is , is not centrally symmetric; whereas [i^K = r~ 2 d/dr is centrally 
symmetric. It is interesting to notice that, apart from a multiplicative constant, 
r~ 2 d/dr is the only incompressible centrally symmetric field in the Euclidean space 
i?3. It is also proved easily that the most general expression for a centrally symmetric 
incompressible field in E n is r x ~ n djdr , where it is understood of course that the 
space has been coordinated with a system of spherical coordinates. 

Another way for obtaining a CRF of the field ( ^I|) is to start with the general 
integral (^4j) and construct an incompressible field L = ^d/dx 1 by formula ([[(]). The 
first component of L is 

f 1 = e tJ V.iV'j = 0,3^,2 - 0,2^,3 = x~ 2 , 
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where = y/x, if) = z/x. By ( |i8| ) \x = ^/f] 1 = x~ 3 is a CRF, which was obtained 
earlier, using a different method. 

The integral curves of K are given by 

dx^x l dt (i = 1, 2, 3). 

Integrating the last system we get x l = x l e*, where x l (t = 0) = x l . The parametric 
equation of the integral curves is r=ro e*. 

The integral curves of the incompressible centrally symmetric field r~ 2 d/dr are 

r = (3t + r 3 y 1/3 } = O , 6 = e . (26) 

Since r > we have t > —Tq/3. When t — > — rg/3, r — > 0, and when i — > +oo, r — > 



+oo. In other words every point moves under the action of the transformation (26) 
in a straight line (0 = 0o, = #o) ensuing from the origin when t —Tq/3 to occupy 
(r , 0o, #o) at t = 0, and goes to infinity as t — > +oo. 

7. A Note on the Algebraic Structure 
We know that the set of all vector fields in M n , named by T, form a Lie algebra 
with respect to addition, multiplication by a scalar and the commutator operation. 
If L = £ l d/dx l and L\ = r/ k d/dx k , then 



L + Li = (C + Tf )d/dx i , cL = (cC)d/dx\ 
[L, Li] = tfrfc-rftyd/to*. 

It is clear that the postulates of a vector space are satisfied and that [ , ] is a commu- 
tator on T. Let 3 C T be the set of all incompressible fields in M n . We shall prove 
that 3 is a sub-algebra of T. Indeed, if L,Li G 3 then 

div(L + Li) = divL + divLi = 0, div(cL) = c divL = 0, 

and hence 3 is a subspace of the vector space T. Furthermore [[! j 

div[L, Li] = L divLi — LidivL, 

and hence div[L, Li] = 0, and 3 is a sub-algebra of T. 

Let II C T be the set of all Killing fields [0] in M n . Since the divergence of every 
Killing field vanishes 0, we have II C If not empty, II is a Lie algebra 0, and 

ncScr. 

Appendix Al. Proof of Theorem (2.1) 
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Necessary conditions are obvious. To prove sufficiency we define a function 

B = — J A\dx 1 dx 2 ....dx n 

Then 

cifi 23 ...n = A 1 , c 2 S 13 ...„ = ( J A 2 2 dx 1 dx 2 ....dx n ), 13 ... n = A 2 (by (6)), etc... 
Appendix A2. Proof of Theorem (3.1) 

By® 

JgdivL = e ll -- l « J B il ....i n _ 2 ,i n _ 1 i n 

For every two values of the indices z„_i and i n there corresponds a unique set of values 
to the indices {h, ..-,1^,-2} such that i\ < .... < i n _ 2 , and hence there corresponds 
to given values of z n _i and i n in the last sum two terms which are equal in absolute 
value: 

^ {Bii i n - 2 ,in-lin ~ -Bi 1 ....i n _ 2 ,i„in-l) ~ 

The multiplier e is equal to +1 or —1 according to the permutation (i 1 ....i n ) is even 
or odd. It follows that the sum over all the values of i n -\ and i n vanishes. 

Appendix A3. Proof of Theorem (3.2) 

We have 

(V9C n ),in = (e il -- <n 01 ) ix^2, ia 0«-l,in-l).*n 

= e ll -- ln (01,i lin 02,i 2 <t>n- X ,i n -x + + 01,ii02,i 2 ^»-l,in-itn) 



lll n 12— ln-2 



+ ( E E ell "' % 01,11 <£»-2,t„J 



Ijl—l^Ti 



We shall prove that the content of every bracket vanishes. For arbitrary fixed values 
of the indices {i^-.-dn-i} the first bracket contains two terms which are equal in 
absolute value and different in sign, and is written as 

e ((j>l,i lin ~ 01 i i n i 1 )(02,i 2 03,i 3 0n-l,i n _i) = 

The other brackets vanish similarly, and the theorem is proved. 
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